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In this paper, we continue the work begun in a previous article. We compute, in the formalism of local com- 
posite operators, the value of the asymmetry in the dimension two condensate for finite temperatures. We find 
a positive value for the asymmetry, which disappears when the temperature is increased. We also compute the 
value of the full dimension two condensate for higher temperatures, and we find that it decreases in abolute 
value, finally disappearing for sufficiently high temperature. We also comment on the temperature dependence 
of the electric and magnetic components of the condensate seperately. We compare our results with the corre- 
sponding lattice date found by Chernodub and Ilgenfritz. 



I. INTRODUCTION 



The dimension 2 gluon condensate (AJJ) in pure Yang-Mills theory has been proposed in [lly], and it has been investigated in 
i 1 different ways since then BI 414I1 . 

In Q3|] an analytical framework for studying this condensate has been developed, based on work carried out in the Gross-Neveu 
model 01511 . Different problems had to be overcome. First of all there is the gauge invariance of this condensate. In order to 
make the operator A 2 gauge invariant, one can take the minimum of its integral over the gauge orbit. Since § d d xA v A u , with 
U € SU(N), is positive, this minimum will always exist. In a general gauge, however, the minimum is a highly nonlocal and 
thus hard to handle expression of the gauge field. A minimum is however reached in the Landau gauge (d^A^ = 0), though, so 
t— I ■ that working in this gauge reduces the operator to a local expression 1 . Secondly adding a source J, coupled to A 1 , makes the 
theory nonrenormalizable at the quantum level. To solve this, a term quadratic in the source must be added, which in turn spoils 
the energy interpretation of the effective action. One way around this is to perform the Legendre inversion, but this is rather 
cumbersome, especially so with a general, spacetime dependent source. One can also use a Hubbard-Stratonovich transform, 
which introduces an auxiliary field (whose interpretation is just the condensate) and eliminates the term quadratic in the source. 
Details can be found in [3]. The result was that the Yang-Mills vacuum favors a finite value for the expectation value of A 2 The 
precise renormalization details of the procedure proposed in yl] were given in Jj]. 

Recently, Chernodub and Ilgenfritz IU2I1 have considered the asymmetry in the dimension two condensate. They performed 
lattice simulations, computing the expectation value of the electric-magnetic asymmetry in Landau gauge, which they defined 

j> ! as 

^ = (8 2 Al)--^- d t(g 2 AJ). (1) 

At zero temperature, this quantity must of course be zero due to Lorentz invariance 2 . Necessarily it cannot diverge as divergences 
at finite T are the same as for T = 0, hence this asymmetry is in principle finite and can be computed without renormalization, 
for all temperatures. 

A first remark concerns the visibility of the (de)confinement phase transition in the value of the asymmetry [12]. At temperatures 
lower than the critical one, the asymmetry goes from zero at zero temperature to a positive value, which reaches a maximum at 
the critical temperature. At higher temperatures, the asymmetry decreases and becomes negative when T > 2.21 T c . The two 
transition points — the phase transition temperature and the symmetric point where the asymmetry goes through zero — divide 
the temperature range in three regions. These seem to coincide with those associated with the condensed, liquid, and gaseous 
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1 We ignore the Gribov problem here, see also [5fl. 

2 We shall deliberately use the term Lorentz invariance, though we shall be working in Euclidean space throughout this paper. 



states of the magnetic monopoles, whose dynamics are closely related to confinement and deconfinement (see, for example, 
IllqO . At yet higher temperatures, one would expect the perturbative behavior to kick in, which goes like 

A A 2(T) = -^8 2 T 2 (2) 

at lowest order. 3 However, lattice artifacts prohibit lattice computations at sufficiently high temperatures to see this 101 . 

At low temperatures, from thermodynamical arguments one would expect an exponential fall-off with the lowest glueball mass 
in the exponent, A ~ e~'"s' r . Instead, the authors of IU2I1 found an exponential with a mass m significantly smaller than m g i. So 
far, there is no explanation for this behavior. 

In 1181 . the authors and collaborators have extended the framework from yj to include the asymmetry A A 2. In that article, 
the potential was computed for T = 0, and no non-trivial value for the asymmetry was found — as is necessary for Lorentz 
invariance. In this paper we extend the computations in order to include finite temperature effects, with the aim of shedding 
more light on the results of 1121 . In section [TI] we give a short review of what was found in [18], which is then continued by a 
computation of the finite temperature effective action in section|III] In section[IV]we find and discuss the minima of the potential, 
the values of the different condensates and their temperature dependence. Section[V]concludes the paper. 



II. PRELIMINARIES 

In IU8I1 the effective action in the presence of a dimension two condensate and of an asymmetry was computed. Since the starting 
point for the calculations when T + are identical, we shortly review the steps taken in 1 18]. 

The starting point to compute the effective potential is the following Lagrangian density: 



r.2 



+ i^ + i^AA + i- W )(AX) , (3) 

where the following couplings have been introduced: 

^1/^161^ 
^ g 2 N \U + 52 16k 2 ) ' W 

- n2 - 1 f l 73 g 2 N \ 
W ~ g 2 N U + 1044 16tiV ' (5) 

to one-loop order. The vacuum expecation values of the a and (p^ v fields are 

(a) = -§<A2>, (6) 

<<fV> = -\U'^- h -fA" K A . (7) 



In order to simplify notations, we set 



m 2 =gd = YJ^gC, (8) 

^v=g<v = 4 rf—[g ( ?t"> ( 9 ) 



3 In 1121 the opposite sign was erroneously found, which seemed to agree with the highest temperatures found in the lattice computations. Given the sign of 
J2), one would expect the qualitative behavior of asymmetry to make yet another turn at higher temperatures. 



which denote an effective mass and an effective mass matrix. With these notations, the condensates as defined with the conven- 
tions of IU2I1 are 
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2 N 

where the Latin index denotes the space components. 

With these givens it is possible to compute the effective action for a space-time independent a' and M^ v using the background 
field formalism. We separate the two fields into a classical part and quantum fluctuations, after which the fluctuations can be 
integrated out. Expanding the resulting path integral over A^ to one-loop order gives 



V eS (o',M^)=L[A ft = 0,o',M l 



N z -l 



ftvl 
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(12) 



where the limit t, — >• is implied, as we work in the Landau gauge. As we are interested in the asymmetry, we parametrize the 
mass matrix as 



n 



M llN =A 
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(13) 
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i.e. we preserve rotational invariance in the spatial part. With this form, the trln in the effective potential can be split into 
different parts, and in the limit ^Owe get 



V e «(a',M^) = L[A /i = 0,a',M l 
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(14) 



III. THE TRACES 

At finite T and in Euclidean space-time, the spectrum of — d 2 is discrete — the eigenvalues are 4-n 2 T 2 n 2 +k 2 where n e Z are 
the Matsubara frequencies and k is the momentum in the spacelike directions. It happens to be convenient to take the second 
and the last terms of (fT4l together (mark that trln —d 2 in dimensional regularization does not vanish for finite T, but it gives a 
constant contribution to the energy), so that we have to compute the following traces: 
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(15a) 
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The sums can be computed using standard techniques. In order to have convergent sums, one first writes (for concreteness, 
consider the first sum): 



52 lnUn 2 T 2 n 2 +k 
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Then the sum can be computed, for example using a formula derived from the product representation for the sine function: 

~ 2 x ~ 1 cothB 1 



ie = en i 



j t 2„2 + 9 2 26 e 2 ' 



(17) 



where we have taken the logarithm of both sides of the first identity, and subsequently taken the derivative with respect to 0. 
Applying this to our sum and performing the indefinite integral in /j, we find 



±T / . A \ \ll 2 + m 2 -^ 

= 21nsinh — 

d-lj 2T 



£ In [4% 2 T 2 n 2 + k 2 + m 2 =21nsinh-i +C , (18) 



where C is a constant of integration. By considering the T — > limit, one can show that it must be equal to 21n2. This gives the 
result that 



;v : 2 A \ „„ f d d - l k , , . , \l k2 + m2 -7TT 
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This can be split into the T = contribution and a finite temperature correction: 



trln ( -d 2 + m 2 - - — - | = trln [ -d 2 + m 2 - ■ 



d-lj V d-\ 
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where we have set d = 4 in the temperature correction. 

The second expression can be computed in an analogous way, except that the numerator in the sum is a fourth-order polynomial, 
and it is thus necessary to split the fraction in partial fractions. We find the temperature correction to the trace to be equal 
to 
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where we have used the short-hand notations a = m 2 + 2k 2 — A/3 and (3 = k 2 (k 2 + m 2 +A 



2T I T^TT I ln I 1 ~exp~ v L - L +ln 1-exp- v L - L | | (21) 



IV. MINIMIZING THE POTENTIAL 

From 1 18] we already know the zero-temperature effective potential. If we add the temperature correction found above, we 
can start the work of searching for minima. As the expressions involved are pretty much unhandleable, we use two strategies: 
expanding in series gives some analytical insight in the low- and high-T behavior, and numerical minimization gives a global 
view of the temperature dependence. For the numerical part, we have used/} 2 = 4.12 A?-—, the value of go' in the non-perturbative 

minimum at zero temperature 13], and N = 2. It is possible to have ft 2 shift as go' gets modified at finite temperature; this, 
however, does not significantly change the results. 



A. Numerical minimization 

Plotting and visually inspecting the potential reveals only one minimum, which coincides with the already known non- 
perturbative minimum at T = 0. One would expect the zero-temperature perturbative solution to become a saddle-point of the 
potential at finite T, but it turns out that this saddle-point can only be found from T = 0.45 A^ onwards. For lower temperatures 
it seems that the saddle-point is located in a region of the parameter space where the effective potential has an imaginary part. 
For slightly higher temperatures, the saddle-point and the non-perturbative minimum merge and from a temperature of 0.67 A-^g 
onwards no solutions to the gap equation can be found anymore. We will say more about this in paragraph llV CI 

The values of the condensates in the non-perturbative minimum are plotted in figureQ] In this figure we have used the sign and 
prefactor conventions of IU2I1 instead of those from |3j], which means that the value of (A 2 ), being the opposite of o, is negative. 
We see that the absolute value of (A 2 ) is slightly lowered by raising the temperature. The asymmetry is positive and rising, just 
as was found on the lattice in fl2(| . Our value for the asymmetry seems to be slightly lower, but as we have only done a one-loop 
calculation, one cannot expect the results to have very high accuracy. 

In figure|2]the values of the electric part and the magnetic part are plotted separately. At T = both are, naturally, equal. When 
increasing the temperature, the electric component goes up, while the magnetic component remains approximately constant. 
This is also what has been found on the lattice] 17]. Similar conclusions for correlations in the gluon condensate were also found 
inlET 
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FIG. 1: The (g Az) condensate (full line) and the asymmetry A A i (dashed line) as functions of the temperature, in units Ajyjg. 
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FIG. 2: The electric (dashed) and magnetic (full line) components of the (Az) condensate as function of the temperature. (The magnetic 
component has been divided by three to be able to compare with the electric component.) 

B. Low temperatures 

Analytically, the limit T ~ can be considered. In order to find the dominant behavior, we proceed as follows: First, the 
exponentials in the integrals are small for small T, meaning that ln(l — expx) ~ — expx. Then, we expand the square roots for 
small \k\. The expansion of the square root will only be valid up to a certain value of \k\, but this is still a legitimate step as the 
greater values of \k\ hardly contribute due to exponential suppression. In the expansion we only keep the terms up to the first 
non-trivial order of \k\, after which the integrals can be easily evaluated. 

We find that the three integrals in the potential have lowest-order behavior 
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respectively. It is clear that, for low T, the second integral will dominate. If we take this to be the first low-temperature correction, 
we find for the asymmetry 
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and (to this order in the temperature) there is no correction to (g Az). If we apply a fit to the low-temperature part of our 
numerical data, the two results are in nice agreement. Mark that, as there is no T 4 correction to (g 2 A 2 ) but there is a positive 
one to A A 2, the magnetic component of the condensate will decrease its value, or increase its absolute value, as can be seen in 
figure|2] This is in opposition to the behavior of the electric component, which only decreases in absolute value. The increase in 
| {g 2 A\) | is small, however, and it is not sure how higher-loop correction will influence this result. 

In IU2I1 it was found that the value of the asymmetry was best described by an exponential form 
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with m = 201(8)MeV. They, however, only had data for T > 0.4 T c . For such higher temperatures, the lowest order in the 
expansion is, of course, not sufficient, and the exponential corrections cannot be ignored anymore. In order to investigate the 



behavior of the asymmetry, the numerical results have to be used again. It turns out to be very difficult to find a fit good enough in 
broader intervals. Given the complexity of the analytical expressions and given the fact that we have only done the calculations 
up to one-loop order, it is not possible to say more about it, however. 

Mark that, given the fact that our model has a mass gap, one would expect an exponential behavior with the effective gluon mass 
in the exponential. Our not finding this is due to the Landau gauge prescription: the last inverse propagator in (Tl4T > does not 
correspond to a simple Yukawa form when the asymmetry becomes nonzero. One should mark that the Landau gauge is singled 
out as the gauge where (A 2 ) reaches its minimum along the gauge orbit, giving this condensate a physical meaning. This is not 
the case with the asymmetry, and as such it is not all that surprising to find a non-exponential behavior for it. 



C. High temperatures 

In order to get more insight in the disappearance of all solutions to the gap equation at higher T, we will expand the effective 
potential in this limit. Herefore it is necessary we return to dT5l > and do the integrations first, receiving, for example, for the first 
expression 

^-1 d-2 „, ,_ lvp +~ / 2 , 2 , , ' ^^ 



2 {A%r -m^-^\l_X %T n +m ^) (25) 

This can be expanded in high T, but one sees that the n = term has to be split off. We find: 
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where C,(s) = Ym=i «? ls trie Ri emann zeta function. The term with i = 2 will give a pole in the d — > 4 limit. The second 
contribution from U5[ can be expanded analogously. 

All together, we find the following high-T expansion: 
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where we have dropped the T 4 term, as it does not depend on the fields in any way and is, thus, irrelevant. Now one has to keep 
in mind that, at high temperatures, one expects the fields to scale with the temperature, and that all terms in the expansion above 
are effectively of the same order in T. However, one expects to have that m 2 ~ A ~ g 2 T 2 , making the above series one in g, with 
the first two terms being of the same order. Solving the gap equation perturbatively yields at lowest order 

r, 13N 1 i N ■, ^ 

m 2 = -—g 2 T 2 , A = --g 2 T 2 , (28) 

and for the condensates: 

, , 9v iV 2 -l , , JV 2 -1 , 7 

(g 2 Al) = ~^-8 2 T 2 , A a2 = —^-8 2 T 2 • (29) 

This is exactly the result one expects from a perturbative computation. 

Going to higher order in this expansion, some subtleties are encountered. First note that the effective gluon masses squared 
are negative, and the next term in the expansion contains square roots of the masses. This gives an imaginary part to the 
potential. Another point of note is the fact that our expansion has effectively become a series in the coupling g instead of one 




FIG. 3: The diagrams giving the Debye mass in hard-thermal-loop resummation. The ghost loop is not necessary 12111 




FIG. 4: More diagrams that need to be resummed: two coming from the a part of the LCO Lagrangian, and again the same two diagrams 
from the (p^ v part of the LCO Lagrangian. The index "£, CO" has been added to remind of the fact that this depicts the two four-gluon vertices 

coming from the LCO formalism. The dotted line is the a or (j)^ propagator. 

in the temperature. This means that, when going to higher order in 1 /T, one has to take into account the effect of higher-loop 
diagrams. These two problems are actually related. 

It has been known for a long time that, at higher temperatures, the perturbation series must be reorganized. 4 In ordinary pure 
Yang-Mills theory, this amounts to giving the timelike gluon a Debye mass m^ — ^g 2 T 2 , which effectively resums the hard 
(high momentum) contributions of the diagrams in figure [3] In our formalism, however, there are four additional vertices. This 
gives rise to four extra diagrams that need to be resummed. They are shown in figure [4] Computing these diagrams, it turns out 
that they exactly cancel the lowest-order contribution from the condensate. 

When doing this resummation, one has to watch out for double counting, which can happen when considering diagrams without 
external lines lf21Tl . However, it turns out that this double counting is put right by the mismatch in symmetry factors in the 
diagrammatic expansion of the vacuum energy. As such, we can proceed without worrying about this. Adding the resummed 
diagrams to the result found in ( f2Tb . we find up to the effective order g 3 : 




(30) 



where, again, terms not containing the fields m 2 and A have been dropped. Again solving this perturbatively, we find: 



, 13 , T 2 m D T \ 3 , T 2 m D T \ 

m 2 = g 2 N[ — + •••). A = — gN\ — + •■■), (31) 

and analogously for (g 2 A 2 ) and A A i. This is exactly what one would expect from perturbation theory. 



V. CONCLUSIONS 

We computed the effective action of SU(A0 Landau gauge Yang-Mills theory in the presence of a dimension two condensate 
and an asymmetry in this condensate. Figure Q] is the main result of this article. We find good qualitative agreement with the 
numerical results of [ 12], with some discrepancies due to different definitions — we define the condensates with the perturbative 
contributions subtracted, whence they vanish in the high temperature, perturbative, regime. The quantitative agreement is less 
excellent, but as our computations are just one-loop and the coupling is not all that small, one may not hope for miracles. A 
two-loop treatment, however, is intractable, even at zero temperature! 22]. 

The low-r behavior seems to be best described by A A 2 = aT 4 , as a naive computation in an Abelian Higgs model would lead us to 
expect[ 12]. The mismatch with the exponential fit found in [ 12] is probably due to their having data only for T > 0.4 T c . At high 
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temperatures it turns out that resummation a la hard-thermal-loop is necessary, and not doing this will give no solutions when 
imposing that the effective action be real. After resumming the necessary diagrams, the perturbative values for the condensates 
are recovered. No non-perturbative solutions are found. The only part of the temperature range where we cannot boast good 
results is around the phase transition. At that point the temperature is already too high to trust a simple one-loop computation, 
and the high-temperature expansion cannot be expected to still yield good results at a temperature so low. 

When this article was in preparation, we learned that lattice computations for the full dimension two condensate and for the 
electric and magnetic components separately at finite temperature have been completed ! 1711 . Qualitative agreement is again 
good. For T < T c , it is indeed found that the electric component shows much more temperature dependence than the magnetic 
component, which is nearly constant in that range. 
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